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In this Letter we present a calculation of the one-particle irreducible density matrix p[x) for the 
one-dimensional (ID) Hubbard model in the infinite U limit. We consider the zero temperature spin 
disordered regime, which is obtained by first taking the limit U ^ oo and then the limit T — > 0. 
Using the determinant representation for p{x) we derive analytical expressions for both large and 
small X at an arbitrary filling factor < £i < 1/2. The large x asymptotics of p{x) is found to 
be remarkably accurate starting from a;sin(27r^i) ~ 1. We find that the one-particle momentum 
distribution function p{k) is a smooth function of k peaked at = 2/cf, thus violating the Luttinger 
theorem. 

PACS numbers; 71.10.-w, 71.27.-|-a 



Recently, we reported results on the one-particle cor- 
relation functions of the continuous ID system of impen- 
etrable spin 1/2 fermions in the spin disordered regime 
It was found that the infrared asymptotic behavior 
of the correlation functions, although consistent with the 
assumption of spin-charge separation, is not adequately 
described by the Luttinger model. This is to be con- 
trasted with the asymptotic behavior of the previously 
studied correlation functions of the infinite U Hubbard 
model in the "antiferromagnetic" ground state, under- 
stood as a limit of the ground state of the Hubbard model 
as J7 ^ oo. In the latter case the Luttinger model gives 
correct predictions H, |^ . 

In this paper we explore the ID Hubbard model 0| 



model (0 in the spin disordered regime was obtained Q ■ 
For the equal time correlation function (O the determi- 
nant representation, given in Ref 7], can be written in 
the following form: 

p{x) = ^ / -F{z)B^^{z) det(/ + V){z). (4) 



Here the function F{z) is 
F{z) = 1 + 
The determinant 



2 - z 2z-l 



(5) 



H = -Y. + h-c.) + E Un^A^^,,, (1) det(/ + V) = dk, . . . dk 



AT 



in the limit [/ — > co. Here ^/'^.o- are fermion fields with 
the spin index a |, and Ux^a = a'^x.a are the local 
fermion number operators. We will concentrate on the 
one particle irreducible density matrix 



X det 



F(fci,fci) ••• V{ki,kN) 



V{kN,ki) ■■■ V(kN,kN)^ 



(6) 



p{x) = (V'ItV'o.t)' 
at an arbitrary fixed filling factor 



(2) 



(3) 



and in the limit T ^ 0. 

The ground state of the model |^ at infinite U is in- 
finitely degenerate with respect to local spin rotations . 
Since the limit T — > is taken after the limit U ^ oo, the 
thermal average () in Q reduces to the average over the 
infinitely degenerate ground state. This is what we call 
an average taken in the zero temperature spin disordered 
regime of the model. 

Recently, the determinant representation for the dy- 
namical correlation functions of the infinite U Hubbard 



is the Fredholm determinant of a linear integral operator 
V with the kernel 

Yd, „^ _ e+(fc)e_(p)-e+(p)e-(fc) 
2tan[i(fc-p)] 

defined on [—K,K] x [—K,K]. Here 

K = 2TTg (8) 

is twice the Fermi momentum. The functions e± entering 
Eq. O are defined as follows 



-ikx/2 



e^(k) = —;=e 
\/tt 



e^{k) = \^e^'^-'\l-z). 
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(9) 
(10) 
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The function B^^{z) is 



and 



K 



dke-{k){i + V)-'^e-{k) (11) 



Consider the contour integral in Eq. Q. According 
to definitions Eqs. Q and 1)10(1 the Fredholm opera- 
tor V is linear in z. This implies that the product 
B--{z) det(/— V){z) is analytic in the complex z-plane. 
Therefore, the integral is given by the residue of the in- 
tegrand at the pole z = 1/2 of the function F{z) 

p{x) = Js-(l/2) det(/ + y)(l/2). (12) 

Consider the short distance behavior of p{x) first. For 
any x the kernel Q can be written as a sum of 2x 
separable kernels (recall that a; is a discrete variable, 
a; = 0,1,2,...) 



V{k,p) 



z - 1 
An 



where 
u„i(k) 



m — 1, 



-*(™-*), m^x+l,...,2x 



(13) 



(14) 



Therefore, det(/ + V) can be expressed in terms of the 
determinant 



det(H-y) =det2:,(I + V) 
of an 2x x 2x matrix V: 



z-1 


( Q 




2tt 




Q 



(15) 



(16) 



Here Q and P are x x x matrices with the entries defined 

by 



Qmn — 

P — 
mn — 

where 



siii[K{m — n)] 
m — n ' 

siTi[K{'m + n — x)] 
m + n — X 



71,171 = I, 
n,m = 1, 



,x (17) 
,x (18) 

(19) 



For B^^{z) one has 

2 sin Kx z — I 



aT(I + V)-ib, (20) 



nx An 

where the 2a;-dimensional vectors a and b are defined by 



2 sin Kn 



2 sin K{n - 2x) 
^/n{n — 2x) 



n = 1, . . . , a; 



, n ~ X + I, . . . ,2x 



(21) 



2 sin K{n — x) 
y/n{n - x) 



. .,2x. 



(22) 



Eqs. H15|) through H22|) combined with H12|) are conve- 
nient for the calculation of p{x) at small enough x. For 
example. 



P(0) 
P(l) 
P(2) 



K 

2^' 
snvK 

~2^' 

sin^ K {2n - K) sin 2K 



47r2 



87r2 



(23) 
(24) 
(25) 



With increasing x the complexity of the exact expression 
for p{x) grows rapidly. 

Next, we calculate the long distance asymptotics of the 
density matrix ^ using the determinant representation 
Technically, the asymptotic analysis will be similar 
to that carried out for the continuous limit of the model 
in Ref. 0. 

To calculate det(/ -I- V) write the difference equation 
for the kernel ITJ: 



V{k,p;x + 1) =e 



V{k,p; 



+ ie-{k; a;)e+(p; x) cos ■ 



P 



(26) 

From this equation it follows that 

det(/ + V){x + l;z) = det(/ + F)(x; z)W{x; z), (27) 
where 

'l + ^B+^(x) ^D^+{x) 
^C+-{x) l + ^A^+{x) 



W{x) = det 



(28) 



and 



Aab 
Bab 
Cab 
Dab 



K 



dkea{k)e-'''{I + V)-'[e'^eb{k)], 



dkea{k){I + Vy^eb{k), 



dkea{k)e'Hi + V)-^eb{k), 



-K 
K 



dkea{k)e-'^{i + V)-^eb{k). 



(29) 
(30) 
(31) 
(32) 



K 



The indices a and 6 run through two values: a, 5 = ±. 

The resolvent operator {I + V)^^ and, therefore, the 
functions (|29|I - H32|I can be found from the solution of 
the corresponding matrix Riemann-Hilbcrt problem |^. 
The scheme of the asymptotic solution of the matrix 
Riemann-Hilbert problem associated with the kernel Q 
is very similar to the one given in It is based on the 
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p{x;K) 



FIG. 1: The constant C{K) in Eqs. igSJ and (gSj plotted as 
a function of K. The plot data were obtained from the com- 
parison of the asymptotic formula 13511 with the numerically 
calculated Fredholm determinant ||SJ. 



non-linear steepest-descend method . The main results 
of the asymptotic analysis are as follows. For z = 1/2 



W{x-z^ 1/2) = 2-^''' 



2 1 



2x 



where 



In 2 



5W{x), (33) 



(34) 



The error term 5W{x) decays as for xsinK > 1. 
Solving Eq. (|77ll with W given by Eq. one gets in 
the large x limit 



det(/ + V)(x) = e'^(^)(sinii') 



(35) 



where C(K) is independent of x. Numerically, exp[C(ii')] 
is close to unity for all K, as can be seen in Fig. ^ The 
exact expression for C(0) is given in Ref. and is, nu- 
merically, equal to 0.0550839 .... This agrees perfectly 
with Fig. n 

The asymptotic formula for the one particle density 
matrix Q reads 



pix) 



r2V2e'='(^)(sini^)'^ 



cosh^ {Kiy/2) 
(2 sin if)- 



— uKx _ - 

e X 



r(-w/2)2 x^+^'' T{ivl2f 



(36) 



with the relative correction of the order of x^^ . The for- 
mula H36|l is the main result of the paper. 

Let us discuss Eq. (|^ . The structure of the corre- 
lation function is essentially the same as for the impen- 
etrable fermion gas [ij, |2, UJl ■ The correlation function 
contains the exponentially decaying factor exp(— t/Xx), 
and factors obeying the power law scaling. The complex- 
valued anomalous exponents do not depend on K or, 
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FIG. 2: Density matrix p{x) plotted as a function of K for 
a; = 1, . . . , 4. The asymptotic result 13(511 (solid line) is in good 
agreement with the exact result (dotted line) even for small 
X sin K. 



equivalently, on the filling factor g. A similar situation 
takes place for the infinite U Hubbard model in the Lut- 
tinger regime: the Luttinger scaling exponents do not 
depend on the filling factor 0. The results for the con- 
tinuous model, impenetrable fermion gas 0, 0, , can 
be recovered by taking the limit ii' — > in Eq. H36() at a 
fixed Kx. 

Formally, the asymptotic formula given in Eq. (|36() is 
valid for xsmK 3> 1. Nevertheless, it is remarkably good 
even for x sin i^T ~ 1 as can be seen from Fig.|21 where the 
exact expression obtained from Eqs. H15|l through H22|) is 
compared with the asymptotics Eq. H36p . 

Finally, consider the momentum distribution function 



p{k)^ e-^^-p{x). 



(37) 



Due to the exponentially decaying term in the asymp- 
totic expression Eq. the function p{k) is continuous 
with all its derivatives for all k. Combining the short 
distance representation Eqs. (|15|I - H22|) and the long dis- 
tance expansion (|36|l we plot p{k) ior Q < k < -k dX 
different filling factors g in Fig. |2| Note that the smooth- 
ness of p{k) in the spin disordered regime, considered 
here, is in contrast with the Luttinger regime considered 
in Ref where dp{k)/dk is singular at the Fermi mo- 
mentum kp = K/2, in accordance with the Luttinger 
theorem [llj. Another peculiarity of the spin disordered 
regime is that dp{k)/dk is peaked around k = 2kF as 
it can be seen in Fig. |31 This can be viewed as a mild 
violation of the Luttinger theorem for this system. 
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Zvonarev's work was supported by the Danish Techni- 
cal Research Council via the Framework Programme on 
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FIG. 3: Momentum distribution function p{k) (thick line) 
and its derivative —dp{k)/dk (dotted line) for different filling 
factors g : (a) Q = 1/8 (h) Q = 1/4 (c) g = 0.45 (d) g = 
0.49. The Fermi-Dirac distribution (thin line) corresponding 
to these filling factors is shown for comparison. The function 
p{k) satisfies p{k) — p{—k). 



Superconductivity. 
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